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Abstract. We report on the analysis of high-resolution optical spectra for 77 subdwarf B (sdB) 
stars from the ESO Supernova la Progenitor Survey. Effective temperature, surface gravity, 
and photospheric helium abundance are determined simultaneously by spectral line profile 
fitting of hydrogen and helium lines, and are found to be in agreement with previous studies 
of sdB stars. 24 objects show spectral signs of a cool companion, being either companion 
absorption lines or a flux contribution at Ha. Five stars with relatively high luminosity show 
peculiar Her profiles, possibly indicating stellar winds. Our results are compared to recent 
theoretical simulations by Han et al. (2003) for the distribution in effective temperature and 
surface gravity, and are found to agree very well with these calculations. Finally we present a 
binary system consisting of two helium-rich hot subdwarfs. 
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1. Introduction 

In the past years, subdwarf B (sdB) stars were subject to many observational 
and also theoretical studies, raising lots of new questions in attempting to 
answer the old ones. How are these stars formed? The relative importance of 
single-star formation as well as binary evolution still needs to be determined 
from observations. Maxted et al. (2001) showed that many sdB stars reside in 
close binaries, indicating a former mass transfer phase to account for the thin 
envelope of the star. Han et al. (2003) examined various formation channels 
for simulating sdB formation in this context, combining them to yield the ob- 
served sdB population. Theoretical studies of this kind need to be compared 
to observational data of high quality to evaluate the simulation results and 
judge our current understanding of sdB formation. With our sdB sample from 
the ESO Supernova la Progenitor Survey (SPY, Napiwotzki et al. 2001), we 
present a homogeneous, high quality dataset of higher resolution and larger 
wavelength coverage than previous studies of comparable size. Its analysis 
should therefore bring new insight into the physics of sdB stars. 

* Based on observations collected at the Paranal Observatory of the European Southern 
Observatory for program No. 165.H-0588(A) and 167.D-0407(A). 

RSj © 2008 Kluwer Academic Publishers. Printed in the Netherlands. 



lisker.tex; 2/82/2008; 7:22; p.l 



2 



2. Observation, reduction, and spectral classification 

Spectra were obtained at ESO VLT with the UV- Visual Echelle Spectrograph 
(UVES) at UT2 (Kueyen). A slit width of 2'/l was used, resulting in a re- 
solving power of 18 500 (0.36 A at Ha) or better. The wavelength coverage is 
3300-6650 A, with gaps at about 4500-4600 A and 5600-5700 A. For most 
of the stars, two spectra in different nights were taken. They were then re- 
duced with a procedure developed by C. Karl using the ESO MIDAS software 
package and parts of the UVES automatic reduction pipeline. The reduced 
spectra were convolved with a Gaussian of 1.0 A FWHM and rebinned to 

o 

0.4 A pixel size. 

After a first selection of potential sdB spectra from SPY by visual in- 
spection, a rough line profile fit was performed, and objects with helium- 
dominated atmospheres were excluded. Following the common assumption 
that sdB and sdOB stars belong to the same evolutionary type, all selected 
stars are named sdB stars, being 77 objects with 20 000 K < T e g < 38 000 K, 
4.8 < log(g) < 6.0 and -4.0 < log(y) < -0.8 (photospheric helium abun- 
dance y = Nyie/Nu)- 51 program stars are objects from the Hamburg/ESO 
survey (Wisotzki et al. 1996), 9 from the Hamburg Quasar Survey (Hagen et 
al. 1995), 14 from McCook & Sion (1999), and 3 from other sources. 

In many spectra, signs of a cool companion can be seen, e.g. a broad 
CaK line (3933 A) as opposed to the frequently occuring narrow interstellar 
CaK line. The Mgl triplet at 5167 A,5173 A, and 5184 A turned out to be a 
good indicator for cool companions, because it reveals even the presence of 
companions with a very low contribution to the total flux. In the spectra of 
five stars, though, only a rise in the continuum level at Ha pointed towards 
a composite object. Altogether, 24 out of 77 sdB stars show signs of a cool 
companion. These objects are excluded from the following analysis, because 
their determined parameters cannot be considered reliable. 



3. Results from spectral analysis 

The determination of T e g, log (g), and log (y) was done by fitting simultane- 
ously each hydrogen and helium line (except Ha) to synthetic model spectra, 
using a procedure by R. Napiwotzki (Napiwotzki et al. 1999) based on Saffer 
et al. (1994). See Edelmann et al. (2003) for a more detailed description. 
The final results were calculated as the mean value of the fit results from 
individual exposures of each star. For six objects, only one useful exposure 
was available. 

The statistical l-tr-errors yielded by the fit procedure are typically lower 
than 100 K, 0.02 dex, and 0.04 dex for T e g, log(g), and log(y), respectively. 
However, since the differences in the fit results of individual exposures of 
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Figure 1. Distribution of our sdB sample in the r cff -log (g)-plane. Grey star symbols denote 
objects with peculiar Ha profiles (see Sect. 3.2). 



each star are mostly larger, we used the distribution of those differences for 
estimating the "true errors", leading to A T e g = 360 K, A log (g) - 0.048 dex, 
A log (y) - 0.045 dex. 

3.1. Atmospheric Parameters 

Figure 1 shows the position of our sdB stars in the r e ff-log (g)-plane, along 
with the Zero-Age EHB and Terminal-Age EHB for solar metallicity (ZA- 
EHB and TAEHB, Dorman et al. 1993) and the Zero- Age Main Sequence for 
pure helium stars (He-ZAMS, Paczyhski 1971). Up to effective temperatures 
of T s g « 33 000 K, most of our stars lie on the so-called EHB strip (defined 
by ZAEHB and TAEHB), while for higher T e g, almost all objects lie near 
the TAEHB or above it. This behaviour is reproduced in current theoretical 
calculations, as discussed in Sect. 4. 

Our distribution of helium abundances (Fig. 2) shows a wide spread from 
log(y) < -4 to slightly supersolar helium abundance. Compared with the 
values from Edelmann et al. (2003) shown in Fig. 2, there is overall agree- 
ment, but we find a higher fraction of stars with very low helium abundances, 
lacking an explanation. Theoretical calculations of the photospheric helium 
content, which include mass loss rates of M = 10" 14 to 10~ 13 M Q /yr, are able 
to explain values of -4 < log (y) < -2 (Fontaine & Chayer 1997; Unglaub 
& Bues 2001). In our data, 18 out of 53 stars (one third) have log (y) > -2, 
clearly pointing towards higher mass loss rates in many cases or alternative 
(additional) physical mechanisms. 
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Figure 2. Shaded histogram: Helium abundance distribution of our sdB stars. Open his- 
togram: Distribution from Edelmann et al. (2003). 
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Figure 3. Two spectra with peculiar Ho- profiles, the left with a flat line core, the right with a 
faint emission. Observed spectra are shown as histograms, model spectra as solid grey lines. 



3.2. Peculiar Halpha profiles 

Five of our sdB stars show abnormal Ha line core profiles (Fig. 3), having a 
lower depth than the model line core, being broader, and in two cases showing 
a faint emission. No significant variation is observed between the first and 
second exposure of the stars (time intervals range from 3 d to 1 yr), but slight 
changes cannot be excluded. All objects are amongst the hottest program stars 
with r e ff > 34 000 K (Fig. 1), and four lie among the highest luminosity 
objects of the sample (the luminosity of the fifth sdB is about average). The 
helium abundances are in the normal range of -3.5 < log (y) < -1.1. 
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Figure 4. Luminosity in Eddington luminosities versus cumulative normalized counts 
from this work as well as from three previous studies. ZAEHB and TAEHB are for 
[Fe/H] = 0.00, -0.47, -1.48, -2.26 from left to right (Dorman et al. 1993). 

Line broadening could possibly be caused by rotation. However, con- 
volution of the synthetic spectra with rotational profiles does not lead to 
acceptable fits for any hydrogen or helium line. Even for Ha itself, the match 
is still not satisfying, implying that other physical effects have to be the cause. 

There is an obvious similarity to the findings of Heber et al. (2003), who 
observed such Ha profiles in four high-luminosity sdB stars and suggested 
stellar winds to be the cause. Their profiles also show broadening around the 
core region, with the core being flattened in one case and containing emission 
in the other cases. As shown by Vink (2003, these proceedings), it is possible 
to reproduce the emission when calculating winds of M - 10" 115 M Q /yr. The 
broadening, however, is not accounted for, and still remains to be explained. 

3.3. Comparison with other samples 

SdB samples of comparable size were presented in the last decade by Saffer 
et al. (1994), Maxted et al. (2001), and Edelmann et al. (2003). They can best 
be compared to our sample by using their respective cumulative luminosity 
distribution, shown in Fig. 4. It can be seen that the SPY data agree well with 
the observations of Maxted et al. (2001) and Edelmann et al. (2003), whereas 
there is a larger average offset to the data of Saffer et al. (1994) (which 
has already been discussed in Edelmann et al. 2003). Different resolution, 
signal-to-noise ratio, homogeneity, and models used for line profile fitting 
can probably explain these moderate differences between the samples. 
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Figure 5. T^-log (g)-distribution of our sdB sample (filled circles) along with the simulated 
stars from the best-fit model of Han et al. (2003, grey dots). 



4. Observation versus Theory 

In order to draw conclusions about potential sdB formation scenarios, we 
now compare our data with the simulations of Han et al. (2003), which aim to 
determine the relative importance of different sdB binary formation channels 
(stable Roche-lobe overflow channel, common envelope ejection channel, 
helium white dwarf merger channel). Several poorly known physical parame- 
ters, e.g. common envelope ejection efficiency, are varied in order to produce 
twelve simulation sets with different parameter configurations. Out of them, 
Han et al. (2003) determined a best-fit model by comparing periods and min- 
imum companion masses for binary sdB stars with observations by Maxted 
et al. (2001) and Morales-Rueda et al. (2003). A comparison with our r e g- 
log (g)-data is interesting because of the homogeneity and high resolution of 
our sdB spectra, surpassing previous observational samples. 

Before comparing theory and observation, one must consider observa- 
tional selection effects, in the present case mainly due to main sequence (MS) 
companions, which either outshine the sdB and thus prevent its detection, or 
make its analysis unreliable. Han et al. (2003) apply this effect by excluding 
all systems with a MS companion having T e s > 4000 K, which is matched 
by our exclusion of all composite spectra and therefore allows us to compare 
our data with the theoretical predictions. 

Figure 5 shows our data along with the simulated sdB stars from the 
best-fit model of Han et al. (2003). The simulations match well the shape 
of our sdB distribution as well as the r e ff-log (g)-area covered. One of our 
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Figure 6. Comparison of first exposure (grey histogram) and second exposure (black his- 
togram) of HE 0301-3039. Both spectra were brought to the same flux level. Barycentric 
corrections have been applied, and the given wavelength is relative to the apparent line center. 



objects is a potential blue Horizontal Branch star (labeled "1" in Fig. 5) and 
is therefore not subject to a comparison with calculations for EHB stars. The 
object labeled "2" lies in a region where there is almost no simulated star. 
However, this object may be in a fast post-EHB phase of evolution (see also 
Fig. 1). As long as theory does not predict the formation of such an object to 
be impossible, it cannot be considered a disagreement with the simulations to 
observe one such star. 

The positions of our stars at or above the TAEHB at higher tempera- 
tures are reproduced in the simulations, partially reflecting their prediction of 
masses in the range 0.4 < M s dB < 0.6 M , as opposed to 0.46 < M S( jb < 
0.50 M Q from single star evolutionary calculations (which define the ZA- 
EHB and TAEHB shown in Fig. 5) . Note, however, that the latter predict 
some 10 % of the objects to be in the post-EHB stage (Dorman et al. 1993), 
which could probably explain the observed stars at high temperatures and 
luminosities above the TAEHB as well. 



5. A subdwarf + subdwarf binary 

Figure 6 shows parts of two spectra of the helium-rich hot subdwarf HE 
0301-3039, the time interval being 10.5 months. The lines split into two 
apparently equal components, being obviously deeper when combined than 
when separated, the radial velocity difference being about 160km/s for all 
lines. Hence HE 0301-3039 is a double-lined binary system consisting of 
two very similar hot subdwarfs. To our knowledge this is the first such system 
discovered. Further observations are scheduled to measure its radial velocity 
curve and system parameters. 
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6. Conclusions 

We have reported on the spectral analysis of 77 sdB stars from SPY, as well as 
on the detection of the first known binary consisting of two hot subdwarfs. 24 
sdB stars show signs of a cool companion. Of the 53 non-composite objects, 
five stars with relatively high luminosity show peculiar Ha profiles, possibly 
indicating stellar winds. Our data are found to be in good agreement with 
previous studies of sdB stars. The best-fit model of the binary population 
synthesis calculations by Han et al. (2003) reproduces the observed sdB dis- 
tribution very well. Note, however, that at this point no statement can be made 
as to the relative importance of single-star formation channels. A detailed 
comparison of our observations with all simulation sets of Han et al. (2003) 
will be presented in a forthcoming paper. 
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